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1. Introduction 

Calculations in chiral perturbation theory [|], ||, || beyond one loop have proven to be 
intricate. One of the difficulties lies in the contributions from genuine two-loop integrals 
involving massive internal lines. In the case of three-flavour chiral perturbation theory these 
can in addition involve more than one mass scale which so far often prohibited the analytic 
evaluation. For example, an analytical expression for the pion mass in this framework is 
not available. In the present paper, we derive an explicit representation for the two-loop 
correction to the pion mass in terms of dilogarithms. Our final result is still not entirely 
analytic because the closed-form expression for the KKr\ intermediate states is missing. 
We account for these contributions in terms of a simple approximate representation suitable 
for all practical purposes including investigations of the quark mass dependence. 

In the two-flavour version of chiral perturbation theory ^] the pion mass to two loop 
order was calculated in Ref. Q and further useful representations of the result can be found 
in Refs. ||] and Q. In this case, the calculation only involves equal mass on-shell sunset 
integrals which were evaluated analytically. The two-loop calculation of the pion mass in 
three-flavour chiral perturbation theory || was performed in Refs. |7|,||. These calculations 
were performed in the isospin limit m u = m^. A calculation of isospin breaking effects in 
the meson masses is described in Ref. Q. For a recent review of chiral perturbation theory 



beyond one loop we refer to Ref. [10] and an extensive list of references can also be found 



in Ref. [11]. A detailed description of the calculational methods used in the framework of 



the two-flavour theory is provided in Ref. ||i~2||. Recently, explicit analytic results involving 
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the strange quark flavour have been obtained for the two-loop matching of the low-energy 
constants B |lT| , F and £%,... ,£7 |l3| . The analogous relation for the strange quark mass 



dependence of the up-quark condensate (0|im|0) was derived in Ref. [14|. In this last case, 
there are no contributions from sunset integrals. 

On the other hand, the final results for the pion mass in Refs. Q, |j involve contribu- 
tions from sunset integrals that have not been evaluated analytically. Among these, there 
are furthermore tensorial sunset integrals of rank 1 and 2. 1 By Lorentz invariance, such 



tensorial integrals may always be related to linear combinations of scalar integrals [15]. 
In the case of a one-loop integral, the scalar products involving the loop momentum in 
the numerator may then always be re-expressed in terms of the denominators, masses and 
scalar products of the external momenta. Thus, the procedure leads to a representation of 
any tensorial integral in terms of the corresponding scalar integral plus simpler functions. 
In the case of the two-loop 2-point sunset integral this procedure does not go through, 
because the number of scalar products that can be formed with the two loop momenta and 
the external one (five) is greater than the number of available denominators (three). In 
general, one thus ends up with additional integrals with irreducible numerators. 



In Refs. |lj| 17] Tarasov has pointed out that such integrals with irreducible numer- 
ators can always be related to integrals with unit numerators, at the price of shifting the 
dimension d in dimensional regularization by multiples of 2. This then calls for finding a 
connection of integrals with unit numerators in different dimensions. In fact, such relations 
can be found easily from the Feynman parameter representation for a given integral. On 
the basis of these observations, Tarasov shows that, up to contributions from tadpoles, any 
(tensorial) sunset integral with denominators raised by arbitrary positive integer powers 
can recursively be expressed as a linear combination of four master integrals, with co- 
efficients formed with the three masses, the momentum and the dimension d. The four 
master integrals consist of the sunset integral with unit powers of the denominators plus 
those three where in each case one of the denominators appears squared. The recurrence 



relations of Ref. [17] have been tested and implemented in a computer program package in 



Ref. [18]. For a review of the 'integration by parts' relations [19, 20, 21] that are at the heart 
of these results we refer the reader to Ref. |2^] which also contains further bibliographi- 
cal references. The problem of (massive) tensorial two-loop integrals is also considered in 



Refs. [E3L 24, 



As an application of the general results of Ref. [17], we derive the reduced expressions 
for the tensorial sunset integrals relevant for chiral perturbation theory calculations. Their 
representation in terms of the four master integrals is provided for the general case with 
arbitrary momentum and three different masses. With these expressions at hand, we then 
work out a simplified representation for the pion mass to two loops in three-flavour chiral 
perturbation theory. Our derivation relies on the representation given in Ref. || where 
the contributions of the coupling constants of the order p e chiral Lagrangian are given in 
the standard basis determined in Ref. [30], for related work see also Ref. [31]. Our explicit 
result shows that the consistent reduction to master integrals leads to a final result in terms 



1 In the two-flavour case, the tensorial integrals were evaluated along with the scalar integral ]l2 
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of fewer basis functions compared to Ref. ||. 

An account of analytically known special cases of two-loop 2-point functions can be 



found in Refs. [32, 33, [34|]. For recent analytic results for the equal mass sunset graph we 
refer to Refs. [35, 36]. For the pion mass we arrive at an explicit representation in terms of 
dilogarithms, involving in the case of the contribution from the KKr] intermediate states 
an expansion for a small pion mass. Our result is presented in such a manner that the 
renormalization scale dependence of the various contributions is manifest: For example 
we separately display the contributions involving double and single chiral logarithms and 
one separate term accounts for the finite and scale invariant contributions from the sunset 
integrals. The advantage of such a representation is exploited in the subsequent numerical 
study of our results, where we in particular investigate the scale variations of the individual 
contributions. Previous numerical investigations of the two-loop corrections to the pion 



mass can be found in Refs. ||, ||, 37, [3^, 39, 4C, 41, 42, 43 1 and a brief overview of these is 
given in Ref. |44]. 

The paper is organized as follows: The recurrence relations and results for the tensorial 
sunset integrals are presented in Section [2]. The dimension d is left arbitrary throughout 
this section. In Section ||, we then provide the relations required to separate the divergent 
and finite parts in the vicinity of d = 4 and give explicit expressions for the finite parts 
for a number of special cases of interest. In Section ^j, these results are applied to yield 
a simplified representation of the pion mass. The algebraic evaluation is followed by a 
discussion of the numerical implications of our results in Section || Finally, Section || 
contains our conclusions. For the purpose of illustration, we provide in Appendix |A| a 
discussion of the recurrence procedure in the case of the one-loop 2-point function. A second 
Appendix |B| demonstrates the absence of spurious singularities in the limit of vanishing 
external momentum in the reduced expressions for the tensorial sunset integrals. 

2. Recurrence relations for the sunset integral 

We consider the family of two-loop integrals defined by 



S d {M 1 ,M 2 ,M 3 ;p 2 } = (2.1) 



i 2 I 



d d k d d l 1 1 



i 2 J (2n) d (2ir) d [Mf -I 2 - i0] a [Mf - (k - I) 2 - iOf [Mf - (p - k) 2 - iQp ' 

These integrals are often referred to as 'sunset' integrals, in allusion to the corresponding 
Feynman graph shown in Figure |l[ Making use of standard techniques, we arrive at the 
following Feynman parameter integral representation for the sunset (see e.g. Ref [|ll]]) 

i r( a + /3 + 7 -d) +,-4-1,0+,- ■ 



SU,,m,M 2 ,M 3 ,,} = j-^ r(a) , (0)m J (2.2) 



r 1 d 

x / dy y a+l} - 2 ~ V -1 [xyMf + xyM% + xxy(M 3 2 - ys) - i0] d ~ a -^ ; 
J o 

where x = 1 — x, etc. The form of this representation suggests the existence of relations 
among integrals with different indices a, (5 and 7. In fact, Tarasov [17| has shown that 
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Figure 1: The sunset topology. 

all the sunset integrals with positive integer indices may be expressed in terms of the 4 
integrals satisfying a + (3 + 7 < 4 plus products of tadpole integrals, 

I d \M\ = l - I ^ 1 = 1 r (»-D M ^ f2 3) 

al ] i J (27r) d [M2 - A;2 - j0] a ( 4vr )f T(q) • 1 ; 

The recurrence relations for the sunset integrals are quite involved and shall not be repeated 
here. For the purpose of illustration we list the analogous relations for the one-loop 2-point 
function in Appendix [A|. For the tadpole, the recurrence relation reads, 

a — - — 1 

ld a { M } = {a _\ )M2 itiiM} , a + \. (2.4) 

For specific configurations of masses and momentum there are additional relations 
leading to a reduced number of master integrals, 

Sl,i,2{M 1 ,M 1 ,M 2 ,M 2 } - -__S lfl>1 {.} -—S 2)lil {-} - 8(d _ 3) M 2 M2 , (2-5) 

Sl hl {M,M,M;M>} = -^=| Sf^-} - J* ~_ |^ if {M} 2 . (2.6) 

Here and in the following, the symbol {•} appearing on the r.h.s. of an equation is a 
shortcut for the repetition of the set of arguments introduced on its l.h.s. The relation 



Eq. fl2.5|) was given in Ref. [EH] and Eq. (|2,6|) was discussed in Ref. [18], but both may also 



be inferred from the general expressions given in Ref. [17|. 

In general, contributions of graphs with indices a, /3, 7 other than unity are not typical 
and chiral perturbation theory is no exception. However, due to the derivative couplings in 
the chiral Lagrangian one commonly encounters tensorial integrals. 2 Below, we will show 
how these may be reduced to the same set of master integrals. It is convenient to consider 
the tensorial integrals associated with the totally symmetric traceless tensors t Pi " '^"(l), 

I 2 

T »(l) = p , T ^(l) = \n v - - , (2.7) 
I 2 

T ^p{i) = m v i p - -—(g^ip + g v n» + g p n v ) , etc. 



2 See for example Refs. 
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For n = 1, 2, 3, ... we then define the functions S^f {Mi,M 2 , M-y,p 2 }, 

r^-^ip) S^ d l {M l ,M 2 ,M 3 ;p 2 } = 

1 f d d k d d l t^-^(1) 1 



(2.8) 



(2ir) d (2ir) d [Mf - I 2 - i0] a [M 2 



I) 2 - iOf [Mf - (p - A;) 2 - i0]7 " 



dHT^-^{l + xp)f{l 2 )=T^-^{ P )x n d d lf(l 2 ) 



The validity of Eq. fl2.Sj ) can be verified, e.g., via the same steps that lead to the Feynman 
parameter integral representation in Eq. Q2.2| ) and by using the relation 

(2.9) 

In doing so, we moreover arrive at a general result for the tensorial integral of rank n, 

(2.10) 



a(n)d r M M M 9-. _ u ,2n r(/3 + n)r( 7 + n) d+2n 



r(/?)r( 7 ) 



This relation exemplifies a general finding of Tarasov jig, 17 1: Tensorial integrals can 
be expressed in terms of scalar integrals with the dimension d shifted by multiples of 2. 
It remains to find a relationship between scalar integrals in shifted dimensions. For the 
tadpoles, these relations are again very simple, 



lf{M} 



M 2 



lf- 2 {M} , d ^ 2 . 



2vr(d- 2) 



(2-11) 



In the case of the sunset, a suitable relation can be obtained from the Feynman parameter 
integral representation in Eq. By inserting the identity 1 = y+(l—x)(l—y)+x(l—y) 



under the integral signs, we deduce [17] 

5^ i7 {M 1 ,M 2 ,M 3 ;s}= (2.12) 
(4vr) 2 [a/3S£ 2 ljm7 {-} + /^ . 

On the basis of this equation evaluated for a + (3 + 7 < 4 and subsequent index reduction 
one can derive the linear relations between the master integrals in d and in d— 2 dimensions. 
Again, the formal expressions given explicitly in Ref. [17| are quite involved. Making use of 
these results and Eq. ( 2.1C| ) we arrive at the comparatively simple relations for the tensorial 
integrals with unit indices, 



S^{M 1 ,M 2 ,M 3 ;s} = 
(d - 2)8 - (d - 3)(2M 1 2 



(2.13) 



Mi 



3(d-2)a 
2 s - Mf + 3(M 2 - Mf 



M 



3(d-2)s 



3(d-2)s 



+ i/fjJil^M,} + J?{M 3 }] - i/f{M 2 }/f{M 3 } , 



? Note that in our notation, the superscript '(n)' does not indicate a momentum derivative. 
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s[ 2) 1 d 1 {M 1 ,M 2 ,M 3 ;s} 



(2.14) 



D 

d? 



{2d(d - l)(d - 2)s 2 + (d - 3)[4(d - l)(d - 4)M X 2 + d(7d - 12)£ 2 , 3 ]s 
- d(d - 3)[4(d - 1)M X 4 - (d + 4)M 1 2 S 2 , 3 + dS| 3 + 12(d - 2)M|M 3 2 ]} fif (ljl {-} 
4DM 2 A M (^-l)^.i + (rf-3)S2,3 ^ { . } 



DM2 



|2[2(d - 2)s + (d + 2)M 2 - (lid - 18)Mf ]A S , 2 

+ dA 2 2 + 3(3d - 4)A 2 3} Sf i2jl {-} + (M 2 

. (5d- 8)S S 1 -dM| - (Id- 12)M% Trfrw , . rw _ , w 
+ £>(d - 2) i ^-f i ^ lf{Mi}lf{M 2 } + (M 2 



M 3 ) 
M 3 ) 



D(d-2) 



(lid 2 - 32d + 24)s + d(5d - 8)M 2 - 2d(2d - 3)S 2>3 , 



ds 2 



lf{M 2 }lf{M 3 } . 



We used the following abbreviations 

d 



D 



6(3d-4)(d- l)(d-2) 



M 2 + M 2 , Aij = M 2 — M 2 



(2.15) 



^s,j — s + Mj , <_». S)J — „ , 

We have checked Eqs. and ( gig ) using the computer program package of Ref. |l8| . 

Note that apart from the increasing length of the resulting expressions, there arise no 
complications in the evaluation of the tensorial integrals of rank n greater than 2. These 
are, however, not required for the envisaged application. In the case of equal masses, 
Mi = M 2 = M 3 = M, the relation in Eq. ( gig ) simplifies to 



Mi 



s[Z{M,M,M;s} = ^Sl 1 A-} 



1 



(2.16) 



This relation was already employed in Ref. j|] in order to achieve a simplification in their 
expressions. 



3. Renormalization and explicit results for the master integrals 

The divergent part of the sunset integral is known for arbitrary values of masses and 
momenta |46|, 47, |48| . For s > 0, it is convenient to define a quantity Sf^jMi, M 2 , M 3 ; s}, 



S?JAM u M 2 ,M 3] s} 



,,2d-8 



(4vr) 4 1 (d - 4) 



2 (M 2 + Ml + M 3 2 ) 



(3.1) 



(d-4) 



-M 2 ln^fln^ 



M 2 f21n^-l) +M|f21n^-l) +M 3 2 (21n 



1 



Mfln^fln^f 



Ml 



1 + 



A/I 



Mi In 2£ ( In ^ - 1 



- In — ^ 

2 fi 2 



As usual, the separation of the divergent parts involves an arbitrary scale /i, which also 
enters in the combination \i x customary in chiral perturbation theory || , 

e- r '(D-i 



4tt 



(3.2) 
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Our definition of the divergent piece in Eq. (3.1) is such that the //-dependence of the 
contributions that explode in the limit d — > 4 is cancelled by a finite term, given in the 
third line, up to contributions of higher order in d — 4, viz. dSf^^^/dfi = 0(d — 4). The 
definition 



S 1 , 1A {M 1 ,M 2 ,M 3 ;s} = lim Sf, - SfjA-} 

d-^A L ' ' ' ' 



(3.3) 



thus implies that the resulting expression Sx,i,i{Mi, M 2 , M3; s} is both finite and inde- 
pendent of the renormalization scale [i. The corresponding finite parts for the remaining 
master integrals are obtained from the above relation by differentiation with respect to the 
masses Mi, M 2 or M3. 

We now turn to a number of special cases for which the analytic representation is also 
known for the finite parts and will be useful in the application below. For two equal masses 
and the momentum on the mass shell of the third mass the result is [49. |50| , 



'1,1,1 



{M, M,m;m 2 } 



(3.4) 



M 2 

W) 



Li 2 (l-r) 



7T 

~6~ 



T 


2 9 " 




In t + — 


, + 2 


4 



lnr- (2 + r) 



7T- 

12 



+ 1 



where Li2(z) denotes the dilogarithm, Li2(z) = Yl^Li ^ (\ z \ — 1) 01- Note that the 
above result does not involve an expansion in the ratio r = m I'M 2 . For such an expansion 
the following representation of the dilogarithm is useful 



2 

Li 2 (l-r) = ^ + £ 



n=l 



lnr 

n 



1 



T n (0 < r < 1) 



(3.5) 



In the case of equal masses, m = M, the result reduces to the simpler expression [p2, 12 



S 1AA {M, M,M;M 2 } 



M 2 

W) 



+ 



For the case of two equal masses and zero momentum, s = 0, one finds 



S w {M,M,m;0} 



M 2 

(47T) 



F(r) + -ln 2 T-(2 + r) 



12 + 2 



as above, r = m 2 /M 2 . The function F(t) is given by 



Fir) = - 



4Li 9 



a -I 
a + l 



1 — o 



I + rr ) + 3 



4 

1 - - 

T 



(3.6) 



(3.7) 



(3.8) 



F(t) is real-valued on the positive real axis. In the range < r < 4 where the variable a 
is purely imaginary it is given by 



F(r) 



Im L12 ( e 



i2 arctan . 



Clo I 2 arctan 



(3.9) 



where we have introduced the Clausen function Ci2(x) (51 
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4. The pion mass to two loops in 3-flavour chiral perturbation theory 

At leading order in the chiral expansion, the masses of tt, K and r\ are given by P] 



1 - 2mB , m 2 K = (m s + rh)B , m 2 = ^{2m s + m)B . (4.1) 

As they stand, the formulas are valid in the isospin limit m u = rrid = rh. The three masses 
satisfy the Gell-Mann-Okubo relation 



2 



™t,= -(Amjt-mi) . (4.2) 
The chiral expansion of the physical pion mass M 2 can be written as 

Ml=ml{l + 5^+5^} + 0{f). (4.3) 

with 5^ = 0{p n ). The corrections of relative order p 2 were determined in Ref. ||, 

<f) = * [mill ~ \m% - 12(n£ + m%)(L\ - 2V % ) - 8m 2 (L\ - 2L\) ] , (4.4) 
where the symbols and denote the chiral logarithms 

r p = -J—ln™i-, P = n,K, V . (4.5) 
z(47r) z fi z 

At the two-loop level, using the representation for the pion mass given in Ref. ||, 4 we 
find the following contributions 

4 4) = 4j [ 5 i°°p + c i°s x i°s + c i°s + c i°g xi, + C U + c LiXL +c Ci ] . (4.6) 



The term is given by the tree contributions of the order p Lagrangian [3C] 

(64F 2 )- 1 c Ci = -m A K [C{, - C r 20 - 3C r 21 ] - m 2 K m 2 n [C r 13 + ±C[ 5 - C[ 6 - 3C r 21 - C r 32 ] (4.7) 
— \m\\2{C\ 2 + C{ 3 ) + C14 + C[ 5 + 3C[ 6 + C[ 7 — 3C[ 9 — 5C 20 — 3C 2 i — 2(C 3 " 1 + C 32 )] , 

while the contributions bilinear in the Lj are collected in the term ci iX Lj, 

c Li x Lj = 64 [{2m 2 K + ml)Ll + m\hl\ [(2m 2 K + m\){L\ - 2L r 6 ) + m\{L\ - 2Lg)] . (4.8) 



Further, the quantity c/^ accounts for the terms linear in the Lj and involving no chiral 
logarithms, 

9(4^) 2 c Li = m% [W4L r 2 + f L 3 + 32L^ + f L r 5 - 64(Lg + L 7 + L r 8 )] (4.9) 
- m 2 K m 2 n [l6L r 2 + f L 3 + 64L^ + f L r 5 - 64(2L£ + 2L 7 + L r 8 )] 
+ ml [36LI + 74LJJ + f L 3 - 40L^ - ±fL r 5 + 80L£ - 64L 7 + A8L r 8 ] , 



4 We have checked that the representation for the pion mass as printed in Ref. [g] coincides with the 
representation given in the preprint version of Ref. jjj], provided the tree contributions from the order p 6 
Lagrangian are identified properly. 

5 Our normalization conventions for the C\ are such that they have mass dimension —2 and the canonical 
large N c behaviour. 
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while the products chiral logarithm times Li are collected in the term q og x l., 

cio g xLi = m K [-32(4L[ + L\ + §L 3 - \L\ - L\ + 4Z/ 6 + 2L r s ) l r K 

- if (4L[ + U 2 + L 3 - ^L\ - LI + 6L r 6 - 3L 7 + L r 8 ) Q 



(4.10) 



- m^m^[16(3LS - 4L r 6 )l r n - %{4L\ + L r 2 + L 3 - f L£ + 3L£ - 6L 7 - 2L r 8 ) I 
+ m*[-8(14Lj + 8L r 2 + 7L 3 - 9L\ - 3L r 5 + 16L£ + 6L r 8 ) l r n 

- |(4L; + Z7 2 + L 3 — L4 + L5 — 2Lg) g . 

The term q og collects the contributions which involve single chiral logarithms 



(4vr) 2 c log = - 7 -^[9f K + Q - ^fp[47Z; + 121 K + f Q - ^[373/; - f^] , 

and ciog x log is given by the sum of the terms bilinear in the chiral logarithms, 

4 

_ ™4 ri75//r \2 59 ;r ;r , 25 f;r\21 , m 7r [337 //r \2 , c/r 7 r 1 nr-\21 
c log X log — m K [-^{Ik) ~ H'-fCS "+" 12 vV J + "oT LlTv'nv + ot vrS ~~ 12 Vr/) J 



2 2 4 

mismz m^ 



+ m|-m 2 [|(C) 2 - + §g + |(2(^) 2 + 5/^/; - f (^) 2 )] . (4.11) 



Making use of the known renormalization scale dependence of the low-energy constants L\ 
and C'l || [54j], one readily checks that the sum of the terms given so far is scale independent. 
Finally, the term ci 00 p accounts for the contributions from the sunset integrals in a scale 
independent manner. It can further be decomposed as follows, 



Cloop 



1 



(47T) 



4 

m K 



97n 2 
432 



23 
96 



, 2 2 

+ m K m 7r 



25tt 2 35 
324 + 432 



41vr k 



4429 



1296 3456 



C-KK-k ~i~ C-r\r\-K CKKr] 



(4.12) 



The indices refer to the particle content in the contributing sunset graphs. Explicitly, we 
have 



CRK-n 



--rfrj-K 



CKKr] 



III 



K 



3m\ 



3m 



8ml 



SK,KA m l) + 



2/4 4 \ 

m K {m K — m*J 



1 



(4vr) 4 

2 



m 



K 



mi 



48 + 8 



+ 



in 



K 



16 



In 



mi 



m 



K 



2m% 
In4 + 1 



2JCA' 



7r(m 2 ) 



m, 
T8 



Sr),ri,ir{fn v ) 



hm A K 



43m 



A" 



8ml 



+ 



+ 



2m^ 
3m?. 



(47T) 



36 
" 54 



17m£ 
72 

, 17 2 2 



SK,K,r,{m^) + 

5mi 



A 



1 



m 



K 



mi 



7T 

144 



1 1 1 

H In 

8 12 



54 



2 

2 m »7 



3m^ 

SK,K,2r,(ml) 



^ m K , 1 2 2 

— — + -m K m n 



m 



A 



5m 



A 



16 



In- 



mi 



m 



A 



(4.13) 



>2A,A. 



17 



In the above equations, we have introduced an abbreviated notation for the sunset master 
integrals 



S a p,bQ,cii(s) = S a ^ c {mp,mQ,m R ;s} . 



(4.14) 
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{i,i,i;i} 


{M, Vp; p} 


{VX, V\, Vp; p} 


{l,l,\/A;p} 


(47r) 4 5 ljl)1 


-4.34240 


-4.03404 


-5.38100 


-4.03912 


(4vr) 4 5 2 ,i,i 


1.32247 


2.16682 


2.19634 


0.935212 


(4vr) 4 5i,i, 2 




-5.49841 


-6.53194 


1.55782 



Table 1: Approximate numerical values for various master type integrals, indices and arguments as 
indicated. The quantity ^fp denotes the ratio of the physical pion and kaon masses, p = M^/M^ ~ 
0.0743454. \/A stands for the corresponding ratio of the r\ and kaon masses, in the approximation 
where the r\ mass is expressed through the Gell-Mann-Okubo relation, A = -|(4 — p) ~ 1.30855. 



Note that the splitting into the various contributions is not unique and neither is the choice 
of the scale independent piece ci oop - Our choice has the property that all the pieces do 
separately have a regular behaviour in the limit ra\ — > 0. For ckk-k and c m7T this follows 



from the explicit result for the sunset integral given in Eq. ( |3.4| ) . For the contribution ckkt\ 
we do not have a closed expression. However, the results from Appendix ^ imply that the 
expansion of this term around ra\ = is determined by the function ^^^{l, 1, -7=; 0} 
evaluated for zero pion mass. Explicitly, we find 



(47T) 



til 



K 



+ 



+ 



CRKri 
43 



43 
144 



In" 



4 95vr 2 425 



+ 



864 
! 35 



+ 



192 
4 



(4.15) 



113 1 

lnp+_F(t) + — ln z - + — In 
72 ' 288 Vs; 24 3 288 3 



17 

144 



In p 



5725 
147456 



F(l) + 



1 

144 



In" 



7T 

54 
4 



4096 3 



15 
" 32 

864 



1217 
9216 



p 2 + 0(p 3 ) 



where /9 = ra\jra\ and the function -F(cc) is defined in Eq. fl3.8|) . On the basis of the 
numerical values given in Table [j] one can verify that the deviation between the full result 
and the approximation is inferior to 10~ 4 for the physical value of p. In fact, even for 
p = 1 the three terms given explicitly still account for more than 99% of the full result. 
We conclude that in the whole range of interest Eq. ( 4.15j ) is also valid for studies of the 
quark mass dependence of ckkt] as needed e.g. for calculations in lattice QCD. 6 

As a check of our result for the two-loop correction given in the present section, we 
have verified that its expansion for rh/m s <C 1 matches the two- loop representation for 
the pion mass in two-flavour chiral perturbation theory [Q, ||, p. Furthermore, we have 
thereby recovered the matching relations for the two-flavour low-energy coupling constants 
B and t\ [|| given in Refs. [11] and [13], respectively. As a final remark we note that in 
the SU(3) limit, m w = m# 
dramatically to yield 



m r , 



mp, the contributions from the sunsets simplify quite 



Qoop 



1457 rap 
384 (4^) 4 



Clog X log 



^mp 2 

7, K 1 P) 



Qog 



1363 rap 
72 (4vr) 2 



lp 



(4.16) 



°We thank L. Lellouch for helpful discussions on the subject. 
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5. Numerical analysis 



We conclude the present paper with a discussion of the numerical implications of our 
results. For an investigation of the size of the various contributions it is advantageous to 
re- write the corrections in terms of the physical masses M n , Mx , M v and the physical pion 
decay constant F w §. We denote the corresponding corrections by the capital letter A 
and write 

Ml = m\{\ + Ai 2 ) + A[ 4 )} + 0(p s ) . (5.1) 

The requirement A^ — 5^ = 0(p 4 ) does not determine the form of A^ uniquely. We 
make the choice 

Af = -L [m 2 l; - ±m 2 l; - i 2 (m 2 + m 2 )(li - 215) - 8M*(li - 2i7 8 )] , (5.2) 

where I7 V and L r n denote the chiral logarithms involving the physical meson masses 

1 M 2 

^ = 2W ln #' P = ^ K ^- (5 ' 3) 



The expression in Eq. ([572 ) is thus is strictly independent of the renormalization scale p,, 

,Aa ? ) = 0. (5.4) 

Note that in the formula above the contribution from the combination L\ — 2Lq is enhanced 
by a factor of about 26 relative to the one from L£ — 2Lg. 

The next-to-next-to leading order correction A^ is given by the contributions from 
5^ plus a shift generated by the difference 5^ — A^ . It allows a decomposition anal- 



ogous to the one in Eq. (4J3). The corresponding numerical 7 contributions are listed in 
Table [2[ To show the variation with the scale we give all values for three different choices 
of the renormalization scale p. For the terms involving the order p 4 couplings constants 
L\, . . . , Lg, we used the numerical values provided in Ref. || and fit 10 of Ref. ||, re- 
spectively. While the two sets of values have been determined on the basis of an order p 
§ and an order p 6 || phenomenological analysis, respectively, the difference between the 
resulting numbers is rather small. The following observations can be made: The tree graph 
contributions bilinear in the Li are small throughout and amount to less than 1.2 • 10 -3 for 
both sets of the Lj and all three scales p. The contribution linear in the Lj and involving 
no logs is given by 

9(4?r 2 )2jP 4 [M 4 K (52Z7 2 + f L 3 ) - MlMl (8L r 2 + §L 3 ) + M* + 37Z7 2 + f L 3 )] . (5.5) 



r In the numerical evaluation, we use the Dashen-corrected values of the physical meson masses, M v = 
M n0 ~ 134.98 MeV, M K = J %{M 2 K+ + M 2 R0 - M 2 + + M*„) ~ 495.03 MeV, M v = 547.51 MeV and F n = 
92.42 MeV 0. 
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LI 


loop 


log x log 


log 


log xLi 


Li 


L{ x L j 


s 


M v 


Ref. § 


7.60 


0.94 


3.33 


-0.56 


2.94 





14.3 


Ref. § 








0.61 


2.93 


0.11 


15.5 


770 MeV 


Ref. § 


7.60 


3.89 


14.0 


-3.91 


0.52 


-0.01 


22.1 




Ref. § 








-0.29 


0.51 


-0.002 


25.7 


lGeV 


Ref. § 


7.60 


8.60 


22.2 


-11.4 


-1.34 


0.07 


25.7 




Ref. § 








-5.88 


-1.35 


0.001 


31.2 



Table 2: Numerical contributions (in units of 10~ 2 ) to the 2-loop correction as defined 
in Eq. (5.1), evaluated for three different choices of the renormalization scale fi. Note that the 
tree contributions from the order p 6 Lagrangian are absent. The other contributions are listed 
separately. For the Li we use the numerical central values taken from Ref. || and fit 10 of Ref. ]{}], 
respectively, as indicated. The last column gives the sum of the contributions in the other columns. 
We repeat that this sum does not account for the contributions of the low energy constants Ci (see 
text). 



For this combination the resulting values of Refs. [|| and || happen to be close despite 
the fact that the values for the individual coupling constants differ considerably. At the 
scale \i = 770 MeV the bulk of the contribution to A^ comes from two-loop contributions 
not involving the L{. Among those terms, the one from the single logs is dominant, and 
the scale independent contributions from the sunset graphs are also sizeable. On the other 
hand, the contribution from the double chiral logs is small - those terms alone only yield a 
poor approximation to the full result. If the two contributions of type log x log and log xLj 
are combined this leads to partial cancellations in both, the size of the contributions and 
their scale dependence. 

We add a remark concerning the contributions from the coupling constants Ci of the 
order p 6 chiral Lagrangian. The contribution to A^ from the Cj consists of three terms 
accompanied by the factors M^, M\M 2 and respectively. The results of Ref. p6| 
imply a non- vanishing contribution only to the third term. In this respect, this contribution 
is similar to the tree result in Eq. (|4.8j ) where the only contribution surviving the large 
N c limit is proportional to L^(L^ — 2L t 8 )/Fq whereas the rest of the terms all involve a 
1/N C suppressed factor of L\/Fq or Lq/Fq. The non- vanishing estimate for the relevant 
combination of Ci from Ref. f56[| involves unknown contributions related to resonance 
matrix elements of scalar and pseudoscalar resonances. On general grounds |57]] one expects 
that the resulting contributions are of the same nature as e.g. the one for (L^) 2 . In view 
of the small ratio M n /Ms, this contribution amounts numerically to very little, 

4M 4 , 

-pP ~ 10~ 3 for M s = 1 GeV . (5.6) 
On the basis of these observations one would not expect a sizeable contribution to A^ from 
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the coupling constants Cj. This is to be confronted with the substantial scale dependence of 
the terms given in Table ^. Assuming that the contributions of the Cj vanish somewhere in 
the range (i = M n ... 1 GeV leads to an estimate of A^ between 0.14 and 0.31, depending 
also on the preferred set of the Lj. 



6. Conclusions 

In the present paper we reviewed the recurrence relations for the sunset integral in the 
general mass case. In Section we in particular provide the explicit representation for the 
tensorial integrals of rank 1 and 2 in terms of the master integrals. In Section ||, we discuss 
the renormalization of the sunset master integrals and provide a collection of known results 
for the finite parts of the master integrals. 

In Section |4], we used these results to obtain a simplified representation of the pion mass 
to two loops in three-flavour chiral perturbation theory on the basis of the result provided 
in Ref. ||]. Our final result involves 6 master integrals accounting for the contributions of 
intermediate KKir, rjrjTr and KKr\ states while the contributions from the intermediate 
3-pion states have been evaluated and are not displayed explicitly. The result of Ref. Q 
was given in terms of 10 functions corresponding to scalar and tensorial sunset integrals. 
Besides the smaller total number of basis functions, our choice has the advantage that a 
nontrivial calculation is only once needed for each intermediate state. Once the unit index 
integrals are known explicitly, the remaining master integrals can be obtained by derivatives 
with respect to the masses. Of those three master integrals the two for the KKtt and r/7/7r 
intermediate states are determined by the same known function given in Section ||. For the 
master integral involving a KKr] intermediate state we do not have an explicit expression. 
However, its expansion for M% — > can be given in terms of the function evaluated for zero 
pion mass (which is also provided in Section ^) and we obtain a numerically very accurate 
approximation of the function. We assume that this representation will be suitable for 
all practical purposes including investigations of the quark mass dependence of the two 
loop contribution to the pion mass. To know the closed-form result would nevertheless be 
interesting. Given our representation this will be the case once the analytic result for the 
sunset integral with two equal masses is known. 

In chiral perturbation theory, the genuine two-loop contributions only represent part 
of the final result at order p 6 . In addition, there arise one loop graphs with insertions of the 
order p 4 coupling constants Li. Further, there are tree contributions bilinear in the Li or 
involving contributions from the coupling constants C, of the order p 6 chiral Lagrangian. 
Unfortunately, the numerical values for the latter are presently not known. In Section [5], 
we show that a simple estimate on the basis of resonance saturation indicates very small 
contributions from these coupling constants. Neglecting them altogether yields an estimate 
of a total order p 6 contribution A^ between 0.14 and 0.31, where the main uncertainty is 
due to the choice of the renormalization scale at which the contributions from the Cj are 
supposed to vanish. Qualitatively, this result for a rather sizeable correction at two-loop 
order is in agreement with earlier findings in Refs. ||, 37, 38, 39, 40, 41, 42, [4^, 44]. 
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Further studies will be required to arrive at a conclusive understanding of the nature 
of the three-flavour chiral expansion. We are convinced that our explicit representation for 
the pion mass at two loop order will be helpful in forthcoming investigations. For instance, 
knowledge of the quark mass dependence of the two-loop corrections will ultimately be 
needed for studies in lattice QCD. Naturally, the methods and results discussed in the 
present paper will also be applicable to other observables in three flavour chiral perturbation 
theory. For a future more comprehensive phenomenological study it would of course be 
desirable to have similarly simplified representations for further of these. 
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A. Recurrence relations for the one-loop 2-point function 



For the purpose of illustration we list in this appendix the recurrence relations for the 
one-loop 2-point function defined by 

1 fj^k 1_ 

1 J (2ir) d [M 2 - k 2 - i0} a [Mf - (p - k) 2 - iOf 
The corresponding Feynman parameter representation is 



li {M 1 ,M 2 ;s}= (A.2) 



<-a,/3 

1 r(a + /3-|) n 



, dxx a - 1 x f3 - 1 [xMl + xMl~xxs-iQ}i-° t - 13 
(4vr)2 r(a)T{/3) J 

with x = 1 — x. In the case of the one- loop integral the recurrence relations follow imme- 
diately from the integration by parts relations associated with the two available momenta 
p and k. Evaluating these for the case of the integral l£ »{Mi, M 2 ; s} in Eq. (AT) leads to 
two equations involving in particular and I d Solving for the former yields 

li +h/3 {M u M 2 ;s} = x{m2 1 m2 g) {(Ml + Ml - s)^^-} (A.3) 
(d-a-2(3)(M 2 -s)-(d-3a)Ml d 20M* d \ 

with A(x, y, z) = x 2 + y 2 + z 2 — 2xy — 2yz — 2zx. The second relation can be obtained by 
the interchange (a, Mi) «-> (/3, M 2 ). The r.h.s of the above equation only involves integrals 
I^ b {-} with a + b = a + 0. Therefore, all the integrals with positive integer indices can be 
reduced to the integral with unit indices and the tadpoles lf{M\} and lf{M 2 }. 
Next, consider the tensorial Feynman integrals I^j^ defined by 



1 f d d k r^-^ik) 



i J (2vr) d [M 2 -k 2 - i0] a [Ml - (p - k) 2 - iO}! 3 



r^"(p)C g d {M 1 ,M 2 ;p 2 }. (A.4) 
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Evaluation via Feynman parameters leads to 

lW{M u M 2 ;s} = (47T)«^^/^{-} • (A.5) 

A relation that allows to decrease the dimension of the integrals can be obtained from 
the Feynman parameter representation in Eq. (A. 2). Inserting 1 = x + (1 — x) under the 
integral sign leads to 

/^{M 1 ,M 2 ; S } = 4vr[a/^ /3 {.} + /3^+ 2 +1 {.}] . (A.6) 

It suffices then to set a = (3 = 1 and to apply the index recurrence relations to the r.h.s. 
to obtain the desired relation 

iffiM^M^s} = J—-L!\(MlMls)I? il {.} + {Ml - Ml + s)!^} (A.7) 



d - 1 8vrs 



(Ml - Ml - s)lf{M 2 } 



Applying the recurrence relations to the tensorial integrals of first and second rank, we 
reproduce the well-know results 

l[f{ Ml ,M 2 ;s} = M? -^ + S lU-} + ^tfW " lf{M 2 }] , (A.8) 

r Wr M M . sT _ d(M? - Ml + S ) 2 - 4M 2 S d d(M* - M* + s) d 

I hl {M 1 ,M 2 ,s}- 4(d _ i)fl2 + 

d(Ml-Ml) + (M-4)s d 

Md^Tj? Il{A/2} - 

Of course, for the one-loop case where the tensorial integrals do not lead to non-reducible 
denominators, there is no need to apply the method described above. We repeat that the 
present appendix is primarily intended to serve as an illustration of the method which in 
the main text is applied to the sunset integrals. It is conceivable, however, that the method 
could be advantageous in an application involving tensorial integrals of high rank. 

Independently of the method used to derive them, the reduced representations for the 
tensorial integrals ([A .8;) appear to involve singularities in the limit s — > 0, while it is clear 



from their definition in Eq. ( A. 4 ) that these must be absent. That this is indeed the case 
can be verified in the following manner: In the case of zero external momentum the master 
integral simplifies to 

T* (M M m ijjM^ - lj{M 2 } 

V {Ml ' M2 ' 0} = Ml - Ml ' (A - 9) 

which is sufficient to prove the absence of the leading singularities in Eqs. ( |A.8|) . The 
coefficients of the non- leading potentially singular terms involve also momentum derivatives 
of the master integral. However, these may again be re-expressed in terms of an integral 
with a shifted dimension 

d n T d SM M . oX , A „\n T ( a + n ) T (P + n ) T d+2n J, /Aim 

cV^ {Ml ' M2 ' s} = (4 " } — ffojr(fl — 7 «+>W-} • (A.10) 

For vanishing momentum, s = 0, they are therefore all expressible in terms of tadpoles and 
one can establish the regular behaviour of any tensorial integral. 
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B. Recurrence relations for the sunset at zero momentum 



If the external momentum vanishes, s = 0, there emerge simplified recurrence relations. 



The relation [17| 



Si +1> sJM 1 ,M 2 ,M 3 ;0} 



(B.l) 



aA(Af?,A*f,Aff) 
x | - [2(a - /3)M| + (d - a - 2/3)(M 2 - M 2 2 - M 2 )] S^ )/J>7 {-} 

— 2/3M| 5 , Q-i i/ 3 + i i7 {-} — ^a^+l.T-li'} 

+ a(M 2 + M| - M|) [S d a+1 ^{-} - ^ +1A7 _ X {-} 

together with its permutations allow one to express any integral with positive integer indices 
in terms of the master integral with unit indices plus tadpoles. Likewise, there is also a 
simplified relation for shifting the dimension jlT]] 



S d + 2 JM 1 ,M 2 ,M 3 ;0} 



(B.2) 



(47r) 2 d(d + 2- a- /3-7) 
A(M 2 , Ml Ml) S d aAl {-} + {-Ml + Ml + Ml) £%_ 1Aj {-} 

+ {Ml - Ml + Ml) S^_ 1)7 {-} + (M 2 + M| - M 2 ) < A7 _!{-} 



Making use of these relations it is straightforward to show that the expressions for the 
tensorial sunset integrals in Eqs. ( 2.13; ) and ( [2.14 ) do indeed have a regular expansion 
around s = 0. Here, the relation for shifting the dimension comes into play because 
momentum derivatives of the integrals may be re-expressed as integrals with a shifted 
dimension: On the basis of the Feynman parameter representation one verifies the following 
general relation for the derivatives with respect to the masses or momentum 







dMl 
(4vr 



77(1 



d 



dMl 



mo 







dMl 



'<<:< / ^ 

ds 



St Aj {M 1 ,M 2 ,M 3 ;s} 



(B.3) 



2n r(a + mi + n)T((3 + m 2 + n)T(j + m 3 + n) ^ d+2n 



r(a)r(/?)r( 7 ) 



a+mi+n,/3+ni2+n,7+m3+n 



{•}• 



Making use of the above results we easily establish the expansion of the master integral 
Sf± ^mK^rriK^m^rn^} around m 2 = 0, 



Si 1A {m K ,mK,m v ,ml} = ^{^{l, 1, ^;0}(1 - ^±|^_J) p ) 



32d 



(B.4) 



3(d-2) 
64d 



If{l} (5d - 9)i?{l} + (d - 3)i?{^} p + 0(p 2 )} 



where p = m^jrr? K . 
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